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Abstract 

The binaries known as cataclysmic variables are particular binary systems in which 
' the primary star (a white dwarf) accretes material from a secondary via Roche- 

lobe mechanism. Usually, these objects have orbital period of a few hours so that a 
detailed temporal analysis can be performed. Here, we present Chandra and XMM- 
Newton observations of a dwarf nova candidate identified in the past by optical 
observations towards the galactic Bulge and labeled as MACHO 104.20906.960. 
After a spectral analysis, we used the Lomb-Scargle technique for the period search 
and evaluated the confidence level using Monte-Carlo simulations. In this case, we 
found that the A-ray source shows a period of 2.03j]o o? nours (3cr error) so that it is 
^ . most likely a system of interacting objects. The modulation of the signal was found 

with a confidence level of >99%. The spectrum can be described by a two thermal 
plasma components with X-ray flux in the 0.3-10 keV energy band of 1.77!tg 19 x 
10 -13 erg s _1 cm -2 . We find that the distance of the source is approximately 1 kpc 
thus corresponding to a luminosity Lx — 2 x 10 31 erg s _1 . 
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1 Introduction 



The observation and subsequent identification of some X-ray sources as low 
mass X-ray binaries (LMXRBs), high mass X-ray binaries (HMXRBs) and 
cataclysmic variables (CVs) plays a crucial role in the study of collapsed ob- 
jects since the physics that they involve is strictly related with the strong 
gravity regime. 



As an example, CVs are binary systems constituted by a white dwarf pri- 
mary star interacting with a secondary object whose spectral type may range 
from G to M classes. The secondary star funnels matter on the primary via 
the formation of a Roche-lobe. Moreover, how the accretion on the white 
dwarf occurs strongly depends on some key parameters as the magnetic field 
strength. In particular, a disk surrounding the white dwarf will form if the 
primary star presents a low magnetic field (as it happens in the dwarf nova 
systems), while a mass stream appears in highly magnetized binaries as the so 
called polars. In the dwarf novae sample, the accreting material may trigger 
some outburst episodes on timescales ranging from weeks to months so that 
the typical brightness of th e system increases of a few magnitudes (see e.g. 



Narayan and Pophamlll993l . an d references therein). A detailed review of the 



CV properties can be found in iKuulkers et al.l (120061 ). 



In the la st decade, different mi cro -lensing programs a s OGLE and MACHO 



see e.g. lUdalski et al.l 119971 . and lAlcock et al.l |2001| ) monitored several di- 



rections in the sky in order to collect data and give proof of the existence of 
dark matter within the Galaxy. These observational campaigns allowed also to 
identify several eruptive variables as CVs and many authors c ompiled detailed 



catal ogues with the properties of the observed sources (see e.g. lCieslinski et al 



2004 for the CVs identified in the MACHO database). 



The cataclysmic variables which are part of the dwarf novae sample are usually 
discovered looking for several observational features in the optical as well as in 
the UV band. However, it has been noted that the efficiency in the detection 
depends strongly on the amplitude of the outbursts (and their separation in 
time) so that observations in bands different than the optical and the ultravi- 
olet ones are required to confirm the nature of these sources. Hence, the X-ray 
emission coming from these source represents a powerful tool of investi gation 



despi te o f the relatively low lu minosity (~ 10 — 10 erg s , see e.g. lLamb 



1982 an dlKuulkers et al.ll2006l ) of dwarf novae. In this respect, XMM-Newton 



(jJansenl l200ll ) t hanks to its l a rge effe ctive area is suitable for this kind of 
search (see e.g. iRamsev et al.l fl2001 al lbh for a detailed analysis of the first 
XMM-Newton observation of a cataclysmic variable system). 



Here, we report on a ~ 100 ks XMM-Newton observation of a dwarf nova 
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candidate towards the galactic Bulge. The target name, as given in the MA- 
CHO database, is MACHO 104.20906.96 and a detailed descr iption of the 
source optical properties can be found in ICieslinski et al.l (120041 ) . The source 
observed by AMM-Newton shows a light curve with a strong periodic behav- 
ior (with a period of ~ 2.03 h). We estimate the A-ray flux of the observed 
dwarf nova to be 1.77 x 10~ 13 erg s -1 cm -2 in the 0.3-10 keV. The paper 
is structured as follows: in Sect. 2, we briefly report about the optical proper- 
ties of the dwarf nova candidate. In Sect. 3 and Sect. 4, we give details on the 
data reduction and on the spectral and timing analysis conducted on MACHO 
104.20906.960, while in Sect. 5 we draw some conclusions aiming that a cross 
correlation search in different bands of this kind of sources is performed in the 
future. 



2 The optical properties of the dwarf nova candidate MACHO 
104.20906.960 



The MACHO project fsee lAlcock et al.lll997l and lAlcock et al.lll999l ) collected 
images of millions of stars in the direction of the galactic Bulge and Small and 
Large Magellanic Clouds to detect micro-lensing events and trace the amount 
of baryonic dark matter distributed in the Galaxy. 

Because of the large amount of data collected in many years, thousand of erup- 
tive variables as well as CVs were identified so that it was possible to compile 
rich catalogues with the main properties of each of the observed source. In 
particu lar, the discovery of 2 8 new dwarf nova candidates was the main re- 

(120041 ). to which we refer for more details about the 



suit of Cieslinski et al. 



characteristics of each of the newly discovered CV systems. 



The result of the analysis of ICieslinski et al.l (120041 ) is listed in their Table 1, 
and the light curves of all the identified dwarf nova candidates are given in 
their Fig. 1. 

In the present work, we study the dwarf nova candidate labeled as MACHO 
104.20906.960 that was serendedipitously observed in 2005. The coordinates 
of the target are a J20 oo = 18 h 05 m 07 s 12 and 5 J2 ooo = -27°43'09"13. In the 
following, we b r iefly report the main results of the analysis conducted by 
Cieslinski et al.l (120041 ) . In particular, they found that the dwarf nova MACHO 
104.20906.960 has a light curve with length 2412 days during which 4 outbursts 
were possibly observed. The source has a typical V magnitude far from the 
outbursts of 17.6 and a median color index (V — R) = 0.69. The magnitude 
at the maximum of the outburst increases of AV = 1.5 (the outburst color 
amplitude being A(V — R) = —0.3). In Fig. [TJ for completeness we give a 
DSS field of view (red band) around MACHO 104.20906.960 (60" x 60"). 
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Here, the red ring is centered on the MACHO database coordinates and has 
a radius of 1" corresponding to the astrometric precision achieved by the 
MACHO experiment. In the same figure we also give a green ring (with radius 
of 2") and yellow ring (with radius of 1") centered on the source identified by 
XMM- Newton and Chandra, respect i vely (for further details see Section 3). 
As clearly stated by Cieslinski et al. ( 20041 ). when the catalogue of the dwarf 
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Fig. 1. The DSS field of view (red band) around MACHO 104.20906.960 (60" x 60") 
is given with the green (2"), red (1") and yellow (1") rings centered on the 
XMM-Newton, MACHO catalogue and Chandra source coordinates, respectively. 

novae candidates was compiled, it was not possible to find a X-ray counterpart 
of any of the newly discovered CVs. In fact, when the authors assumed a 
F x to F opt ratio lower than 0.5, they expected to have all the sources below 
the detection limit (F . 5 _ 2 .5 k ev - 10" 13 erg s" 1 cm ; 2 ) of the ROSAT/P SPC 



observations. Hence, as the same authors suggested (jCieslinski et al.ll2004l ). X- 
ray observations of dwarf novae necessarily require more sensitive instruments 
as the CHANDRA and XMM-Newton satellite. 



For this reason, in the next Section, we will consider a deep XMM-Newton ob- 
servation towards the identified dwarf nova candidate MACHO 104.20906.960 
and use the Chandra sate llite and its better astr ometric precision to confirm 
its nature as predicted by I Cieslinski et al.l (120041 ) . 



XMM-Newton and Chandra observation of the dwarf nova can- 
didate MACHO 104.20906.960 



Here, we are reporting about a XMM observation conducted in the direction 
of the Galactic Bulge and the serendipitous discovery of the X-ray counter- 
part of the dwarf nova candidate MACHO 104.20906.960. The field of view 
around the source was observed by XMM-Newton in two occasions: in Septem- 
ber 2005 (Observation ID 30480) with both the EPIC MOS and PN cameras 
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Fig. 2. A combined image (in the band 0.2-10 keV) of the MOS 1, MOS 2 and PN 
cameras of the full XMM-Newton field of view (left panel) is shown. On the right 
panel, the dwarf nova candidate is indicated by the red circle (with radius ~ 20"). 



([Turner et al. Il200ll ; IStriider et al. Il200ll ) operating in the medium filter mode 



(exposure time of ~ 10 ks) and in October 2005 (Observation ID 30597) with 
the thin filter and for a scheduled observational time of ~ 100 ks. Unfortu- 
nately, in the first observation, the source is close to the camera borders so 
that, in the following, we do not consider it in our analysis. 

The EPIC data files (ODFs) were processed using the XMM-Science Analysis 
System (SAS version 7.0.0). The raw data were processed using the latest 
available calibration constituent files. We ran the emchain and epchain tools 
and obtained several event list files on which the usual screening criteria were 
applied. Thus, after we searched and rejected periods affected by soft protons 
flares the remaining good time intervals resulted in exposures with effective 
time of ~ 96 ks, ~ 98 ks, and ~ 94 ks for MOS 1, MOS 2, and PN, respectively. 



By using the edetect-chain tool, the coordinates of the source MACHO 104.20906.960 
were determined to be a J20 oo = 18 h 05 m 07 s 20 and 5 J20 oo = -27°43'08"76. In 
this respect, note that a systematic shift between the X-ray and optical po- 
sition of the sources of —2.52" and —3.09" (at a level of la) exists in right 
ascension for MOS 1 and MOS 2. The shift on declination in the same two 
cameras is —1.19" and 0.41", respectively. For this reason, in the following we 
assume that the error associated to the two celestial coordinates is at least 

nil 
~ I . 

In Fig. [21 we give a combined image (in the band 0.2-10 keV) of the MOS 1, 
MOS 2 and PN cameras of the full XMM-Newton field of view (left panel). 
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The deep exposure has shown the existence of several new X-ray sourceJ"*"] 
(as the bright source at the up-right corner which was recently identified as 



a M- type X-ray binary with a 494 s-pulse period neutron star, iNucita et al. 



20071 ) which will be investigated elsewhere. On the right panel, the dwarf nova 



candidate studied in the present paper is indicated by the red circle (~ 20") 



The Chandra satellite has a better angular resolution (~ 1", see lThe Chandra X-ray center website 
for details) than the XMM-Newton telescope and therefore it is worth using 
Chandra to improve the X-ray astrometry and remove any doubt that the 
source studied here is the X-ray counterpart of the dwarf nova found in the 
MACHO catalogue. The Chandra satellite observed toward the target in 2003 
(Observation ID 3789) with the ACIS-S camera and for ~ 10 ks . We created 
images in the full 0.3-8 keV band and search for discrete sources by the CIAO 
celldetect tool with a threshold signal-to-noise detection value of 3. This re- 
sulted in the identification of a source (with signal-to-noise ratio of ~ 4) with 
coordinates a J20 oo = 18 h 05 m 07 s 067 and 5 J20 oo = -27°43'08"90 and astrometric 
precision of ~ 1". 

The green and yellow rings in Fig. [1] are centered on the XMM-Newton and 
Chandra source coordinates and have a radius corresponding to the astromet- 
ric precision of each satellite, respectively. Hence, we believe that the XMM- 
Newton identified source corresponds to the X-ray counterpart of the dw arf 



nova candidate MACHO 104.20906.960 detected bv ICieslinski et all fl2004h . 



3. 1 Spectral analysis of MA CEO 1 04 . 20906. 960 



The source spectra were extracted in a circular region centered on the nominal 
position of the target in the three EPIC cameras (extraction circles with radius 
~ 30"), while the background spectra were accumulated in circular regions on 
the same chip and at the same vertical location. The resulting spectra were 
rebinned to have at least 25 counts per energy bin. 

The spectra were simultaneously fitted with XSPEC (version 12.4.0). A single 
temperature thermal plasma model with absorption by neutral gas (MEKAL 
and WABS models in XSPEC) gave a fit with xl = 1-2 (144 d.o.f.) with 
temperature kT = 6.43jli;i7 keV and a hydrogen column density Nh = (1.6 ± 
0.3) x 10 21 cm -2 . The normalization of this component is N = (8.3 ± 0.4) x 
10~ 5 . Adding a second thermal plasma model yielded a significantly better 
fit (xl = 1-00, for 142 d.o.f.) with values of the temperature kT x = 1.27l°:32 



1 This XMM-Newton field of view was also particular interesting because of 
the discovery of the long parallax microlensing event MACHO-96-BLG-5 (see 



Poindexter et al 



Nucita et al 



2005 for a detailed description of the event parameters, and 



20061 for an X-ray search for the lens). 



6 




Energy (keV) 



Fig. 3. The best-fitting model (see text for details) is given together with the MOS 
1, MOS 2 and PN data. 



keV and kT 2 = 10.65lf3o keV. The fitted hydro gen column densitj{EI is Nh = 
(1.81qJ) x 10 21 cm -2 . The normalizations corresponding to the two components 
are N x = (9.1+f J) x 10~ 6 and N 2 = (7.5±gJ) x 10" 5 , respectively. The addition 
of a third thermal plasma component or a blackbody one did not improve the 
fit in a significant way. In Fig. [3J we show the MOS 1, MOS 2, and PN spectra 
(0.3 — 10 keV) for the source and the respective best fits. The absorbed flux in 
the 0.3-10 keV is F^^io = l-48jTo;i7 x 10~ 13 erg s _1 cm~ 2 corresponding to a 
corrected (unabsorbed) flux (in the same band) of Fq^_ 10 = 1.77^o Jg x 10~ 13 
erg s" 1 cm -2 (all the errors are quoted at the 90% confidence level) 

The source count rate detected by the Chandra satellite during the 2003 ob- 
servation corresponded to ~ (4.04 ± 0.02) x 10~ 3 counts/sec. In spite of the 
fact that the count rate is slightly low, we extracted the spectrum from a 
circular aperture centered on the source and, once imported within XSPEC, 
we fitted it by using the same best fit model obtained with the XMM- Newton 
data. After extrapolating the flux in the 0.3 — 10 keV band we got that the 
Chandra absorbed flux is -Po3-iokev ^ 2.9 x 10 -13 erg s _1 cm -2 , i.e. consistent 
with the result obtained by using XA/fM-Newton and quoted above. 



Note that the column density obtained from the fit procedure is consistently lower 



than that estimated to be in our Galaxy (Njj 
Il990h . 



3.3xl0 21 cm 2 . lDickev Lockman 
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3. 2 Timing analysis of MA CEO 1 04 . 20906. 960 



We extracted the light curves of the source from all the three EPIC cameras 
using apertures ~ 30" in radius and centered on the MACHO 104.20906.960 
position. The background light curves were extracted from the same CCD 
where the source was detected, scaled and subtracted from the source light 
curves. To increase the signal to noise ratio, the counts from the MOS 1, MOS 
2 and PN cameras were summed. We repeated the entire process for the energy 
range 0.2-10 keV (to search for any periodicity in the full light curve) and for 
the ranges S = 0.2 — 1 keV and H = 1 — 10 keV to evaluate the hardness ratio 
light curve in the two above energy bands and search for any change in the 
source state (see next). 



The 0.2-10 keV light curve of MACHO 104.20906.960 during the 100 ks ob- 
servation is given in Fig. H] as binned at 700 s. We searched for periodicity 
between 20 sec and 1 h, and found a periodic s ignal at 2.03 h using a Lomb- 
Scargle periodogram (jLomblll976l ; IScargld Il982l ) . The confidence level of our 
detection was estimated with Monte Carlo simulations assuming a null hy- 
pothesis of white noise, and the results are plotted in Fig. with the 68%, 
90%, and 99% confidence levels represented by the full, dotted and dashed 
horizontal lines, respectively. W e found that the 2.0 3 h period is significant 
at a confidence level > 99% (see ICarpano et al.l 120071 for more details on the 



applied technique). Thus, at 3a confidence level, the detected period of the 

h. Assuming that the X-ray emiss ion comes from 



observed source is 2.03i[J;o7 



the b oundary layer (where the X-ray emission is produced, see lKuulkers et al. 



20061 for a review), the period estimated above is consistent with the typical 



dwarf nova orbital period. The XMM light curve folded at the detected period 
of 2.03 h is shown in Fig. [6J 



T 




5 10 15 20 25 

t(hr) 

Fig. 4. The 0.2-10 keV light curve of the 100 ks observation of MACHO 
104.20906.960, as binned at 700 s. Here the MOS 1, MOS 2 and PN counts are 
summed and background corrected. 
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Fig. 5. Search for periodicities in the XMM-Newton light curve of MACHO 
104.20906.960 using a Lomb-Scargle periodogram analysis. The full, dotted, and 
dashed lines represent the 68%, 90%, and 99% confidence levels, respectively. 



Fig. 6. The 0.2-10 keV light curve folded at 2.03 h using 10 bins. Phase zero is 
associated to the beginning of the XMM-Newton observation. 



Using the light curve extracted in the soft (S) and hard (H) parts of the 
spectrum (as described above), we defined the usual hardness ratio HR = 
(H — S)/ (H + 5*) and built the corresponding hardness ratio light curve to 
search for variability in the state of the source. The hardness ratio light curve 
is given in Fig. [7J Here, the solid horizontal line represents the value of the 
average hardness ratio, while the almost superimposed dashed line is the re- 
sult of a linear fit to the binned data, clearly showing that the state of the 
source remains constant (i.e. without overall significant changes) during the 
observation, as it would be expected since we observed the source during its 
quiescent state (i.e. far from an outburst). 



9 



1.0 - 



5 10 15 20 25 

t (hr) 

Fig. 7. The hardness ratio HR = (H - S)/(H + S) light curve (binned at 1800 s) of 
MACHO 104.20906.960. The used energy bands are S = 0.2 - 1 keV and H = 1 — 10 
keV, respectively. 

4 Discussion and conclusions 



In the present work, we studied the X-ray counterpart of the dwarf nova 
candidate labeled as MACHO 104.20906.960. 

A two temperature thermal plasma model with absorption by neutral gas 
(MEKAL and WABS models in XSPEC) gave an acceptable fit with values 
of the temperature kT x = 1.27lg;^ keV and kT 2 = 10.65tf keV. The fitted 
hydrogen column density is Nh = (1.8 ± 0.3) x 10 21 cm -2 cm -2 and the 
corrected (unabsorbed) flux resulted to be Fq^ 10 = 1.77io;ig x 10~ 13 erg cm -2 
s -1 , all the errors being at the 90% confidence level. We found a periodicity 
of 2.03^07 h in the source light curve by using a Lomb-Scargle periodogram. 

If the X-ray measured hydrogen column density (Nh — 1.8 x 10 21 cm -2 ) 
is entirely due to interstellar extin ction, then using the standard Nh — A v 



relation (see e.g. iGorensteinl Il975l ). one has A v ~ 0.9, being this un upper 
limit to the extinction since many dwarf novae show evidence for intrinsic X- 
ray absorbtion. Hence, we estimate the unabsorbed V magnitude of the source 
to be ~ 16.7 

One of the typical indicator used in the study of CV systems is the ratio 
between the X-ray and optical flux, i.e. I = Fx/Fy. In this case, from the 
above estimated V magnitude the optical flux can be readily obtained to be 
Fy — 7.6 x 10~ 13 erg cm -2 s _1 and, by using the measured X-ray corrected 
flux, one obtains / ~ 0.23 which is a typical value for a non-magnetic CV (for 
which it is expected to have I in the range 10~ 3 — 10). This is consistent with 
the hypothesis that MACHO 104.20906.960 is a dw arf nova (non-magnetic) 
system (as first pointed out by Cieslinski et al. 20041 ). 
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Furthermore, it is well known that for dwarf novae a relation between the 
absolute magnitude in the V band at the outburst and the orbital period 
exists (see 



Warner! Il995l for details), i.e. 



burst 
V 



5.74 - 0.259P(h) for P(h)<15h. 



(1) 



In the case of MACHO 104.20906.960, the dereddend burst apparent magni- 
tude (my urst ) and the outburst absolute magnitude (My urst , as given by the 
previous equation once P{h) is known) read out to be ~ 15.2 and ~ 5.21 
respectively, so that the distance to the source can be easily determined to 
be approximately d ~ 993 pc. Hence, assuming isotropic emission, the X- 
ray luminosity of the source (in the 0.3 — 10 keV band) can be derived as 
Lx = 4: jrd 2 Fx obtain ing Ly ~2x 10 31 erg s" 1 , again typical for a dwarf nova 



se e.g. IWarnerlll995l and lKuulkers et al.ll2006l ). 



We e mphasize that long term campaig ns, as the past OGLE (jUdalski et al. 
19971 ) and MACHO (lAlcock et al.ll200ll ) experiments, are particularly attrac- 
tive due to their capability of monitoring several million of stars for long time. 
These observational campaigns, among the other results, allowed the identi- 
fication of several variable sources as CVs so that more studies based on an 
enhanced statistic may be performed. Of course, thanks to the nature of many 
of these sources (interacting binary systems), one can get more information 
if a multi-wavelength search (spanning from the optical to the X-ray band) 
is done as in the case of MACHO 104.20906.960 whose nature is confirmed 
by XMM-Newton data. We suggest that a cross correlation search in different 
bands for this kind of sources is performed in the future. 
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